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ABSTRACT 



Context. Monitoring of pulse period variations in accreting binary pulsars is an important tool to study the interaction between the 
magnetosphere of the neutron star and the accretion disk. While the X-ray flux of the brightest X-ray pulsars have been successfully 
monitored over many years (e.g. with RXTE/ASM, CGRO/BATSE, Swift/BAT), the possibility to monitor their pulse timing proper- 
ties continuously has so far been very limited. 

Aims. In our work we show that the Swift/BAT observations can be used to monitor coherent pulsations of bright X-ray sources and 
use the Swift archival data to study one of the most enigmatic X-ray pulsars, Hercules X-1. A quasi-continuous monitoring of the 
pulse period and the pulse period derivative of an X-ray pulsar, here Her X-1, is achieved over a long time (> 4yrs). We compare our 
observational results with predictions of accretion theory and use them to test different aspects of the physical model of the system. 
Methods. In our analysis we use the data accumulated with Swift/BAT starting from the beginning of 2005 (shortly after launch) until 
the present time. To search for pulsations and for their subsequent analysis we used the count rate measured by the BAT detector in 
the entire field of view. 

Results. The slope of the correlation between the locally determined spin-up rate and the X-ray luminosity is measured for Her X-1 
and found to be in agreement with predictions of basic accretion torque theory. The observed behaviour of the pulse period together 
with the previously measured secular decrease of the system's orbital period is discussed in the frame of a model assuming ejection 
of matter close to the inner boundary of the accretion disk. 
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1. Introduction 

The persistent accreting pulsar Hercules X-1 was one of the 
first X-ray sources discovered by the Uhuru satelhte in 1972 
jTananbaum et al.lll972l : iGiacconi et al. I [1971 and since then it 
remains one of the most intensively studied X-ray pulsars. The 
basic phenomenological picture of Her X-1 was established soon 
after its discovery: a close binary system consisting of an ac- 
creting magnetized neutron star with a 1 .24 s s pin period a nd a 
stellar companion HZ Her (first suggested by Liller 19 72|) - a 
main sequence star of the spectral type A/F ( Crampton 1 1 9741) . 
The mass of the optical companion is ~2Mq which places the 
system in the middle between high and low mass X-ray binaries. 
Other main parameters of the binary system are the following: 
orbital period f oit -1-7 days. X-ray luminosit y of the source 
Lx ~ 2 X lO^^ergs ' for a distance of ~7 kpc ( Reynolds et all 

2009bl) 



Il997h . The binary orbit is almost circ ular (IStaubert et al. 
and has an inclination / ~ 85 - 88° dGerend & Bovntor 



wm . 



The magnetic field strength on the surface of the neutron star is 
believed to be around 3 x 10'^ G, as estim ated from the energy of 
the cyclotron resonant scattering feature dTriimper et alJ[T978l) . 

Like many other X-ray pulsars Her X-1 exhibits signifi- 
cant variation of the pulsation period (i.e. spin rate of the neu- 
tron star). Alternation of spin-up and spin-down episodes on 
time-scales of several months in this syst em is superimposed 
on a background of systematic spin-up dShefi^er et al.l 119921: 
IStaubert et al. 2006; Klochkov 2007). The behaviour of the pul- 
sar's spin period on shorter time scales is not very well stud- 
ied because such a study would require a continuous monitoring 



of Her X- 1 with a sensitive X-ray detector Only between 1991 
and 2000 the BATSE instrument onboard CGRO (iFishman et al.l 
il989 ) provided the information about t he source's pu l se per iod 
on a regular basis. These data allowed [Staubert et al.l (l2006l) to 
reveal an anticorrelation between the pulse period and times of 
X-ray tum-ons, i. e. switching from an off-state with low X-ray 
flux to the so-called main-on state with high flux (such turn-ons 
regularly occur in Her X-1 with a period of ~35 days and are be- 
lieved to be caused by a precessing tilted accret ion disk around 
the neutron star, see e.g. lGerend & Bovntonlll97 6). 

In this work we present a continuous monitoring of the 
Her X-1 pulse period P and its local (measured at the time 
of the observation) time derivative P using the Swift/BAT in- 
strument starting from 2005 (begin of scientific operation) to 
2009. We compare the observed pulse period development with 
that measured previously with CGRO/BATSE. The data of the 
monitoring allowed us to explore for the first time the correla- 
tion between the locally measured P and the X-ray flux of the 
pulsar and compare the results with predictions of the accretion 
theory. The observed strong spin-down episodes are discussed in 
the frame of a model assuming ejection of matter from the inner 
part of the accretion disk along the open magnetic field lines. 

2. Observations 

For our analysis we used the public archival data obtaine d with 
the B urst Alert Telescope (BAT, 150-150keV, Barthelmy et all 
|2005|) onboard the Swift observatory (Geh rels et al.ll2 004). With 
its large field of view ( 1 .4 sterad) the BAT instrument is origi- 
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nally designed to provide fast triggers for gamma-ray bursts and 
their accurate positions in the sky (~4 arcsec). Following such a 
trigger, the observatory points in the direction of the burst, which 
can be then observed with the X-ray and UV/optical telescopes 
onboard the satellite. While searching for bursts, BAT points at 
different locations of the sky, thus, performing an all-sky moni- 
toring in hard X-rays (measurements of the X-ray flux are pro- 
vided by the Swift/BAT team in the form of X-ray light curves 
for the several hundred bright persistent and transient source^. 

Most of the BAT observations are stored in the form of de- 
tector plane maps (histograms) accumulated over the 5-minutes 
exposure time which limits the possibilities of timing analysis. 
In addition to the detector plane histograms, however, the stored 
data contain 64-msec count rates corresponding to the total flux 
detected by BAT. If a bright pulsating source with a known pe- 
riod falls into the field of view of the instrument, the total count 
rate can be used to search for coherent pulsations of that source. 
We have used this strategy to measure the 1.24 s pulsations of 
Her X-1 during its main-on states (when the X-ray flux of the 
source is high) repeating every ~35 days. 




0.5 1 1.5 2 



Phase 

Fig. 1. An example of a Swift/BAT pulse profile used to study 
pulse shifts for the precise measurement of the pulse period vari- 
ations in Her X-1. The profile is taken at iVIJD 54516.083 with 
~ 1200 s integration time. 



3. Timing analysis of BAT data 

To determine the pulse period, we us ed the method similar to that 
described by IStaubert et al.l (l2009bl) . The method includes two 
techniques for the determination of t he period: epoch folding 
with s earch (e.g.lLeahv et alJll983l) and pulse phase connec- 
tion (e.g. lDeeter et al.ll981 ). The first one is used to establish the 
presence of the periodic signal from Her X-1 in the BAT data, 
determine the approximate period, and construct pulse profiles 
(by folding the data with the found period), while the second is 
subsequently applied to the pulse profiles to determine the pre- 
cise value of the period and its time derivative. 

As already mentioned, for our analysis we used the total 
count rates measured by BAT with a time resolution of 64 msec. 
All times of the count rates were translated to the solar system 
bary center and corrected fo r binary motioii (using our newest or- 
bital ephemeris presented in lStaubert et al.l2009 b). Then we per- 
formed a period search using epoch folding in a narrow period 
interval around the expected pulse period (~ 1.237 s). If a strong 
periodic signal is present we determined the period and used it 
to construct X-ray pulse profiles for subsequent puJse phase con- 
nection. The integration time in each case was chosen to be large 
enough to obtain a pulse profile with sufficient statistics (nor- 
mafly ~1 ksec). A typical profile used in our phase-connection 
analysis is shown in Fig.[T] 

To obtain pulse arrival times for the subsequent phase con- 
nection analysis we have fit all individual pulse profiles with a 
template profile constructed by superposing the individual pro- 
files from the same main-on state. In many cases we ignored 
the data at the start and the end of a main-on where the pro- 
files deviate noticeably from those obtained in the middle of the 
main-on. Thus, within one main-on we restricted our analysis 
to the time intervals where the shape of the pulse profile does 
not change appreciably. Nevertheless, we cannot completely ex- 
clude possible systematics which might affect our results due the 
variation of the profile shape (see also Discussion). If the inter- 
val between the adjacent pulse arrival times is short enough (or if 
the uncertainty of the assumed pulse period is sufficiently small), 
one can reconstruct the number of pulsation cycles passed in be- 
tween excluding any mis-counting. In this case, the estimated 
pulse arrival times can be analytically modelled. For example, if 
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Fig. 2. Estimated (observed) minus calculated (assuming a con- 
stant period) pulse aiTival times of Her X-1 in units of its pulse 
period during one of its main-on states observed with Swift/BAT. 
Parabolic fit to the data shown with the solid line corresponds to 
a constant positive P ^ 1.4 x 10"'^ s/s (according to Eq.|2l). 

the pulse period P is constant, the expected arrival time of pulse 
number n is 

tn^to + nP, (1) 

where fo is the arrival time of the "zero"-th pulse. In case of non- 
zero first and second deriv atives of the pulse per iod, the arrival 
times are given by (see e.g. IStaubert et al.ll2009bl) 

t„ = fo + nPo + \n^PP() + \n^PlPo + (2) 
Z 6 

where Pq, Pq, and Pq are the pulse period and its time derivatives 
at the time to respectively. 

A convenient way to explore the variation of the pulse pe- 
riod, often used in phase connection technique, is to construct 
the so-called (O - C) diagram showing the estimated (observed) 
pulse arrival time minus the calculated one assuming a constant 
period (i.e. using Eq.[T]). An example of such a diagram mea- 
sured with JBAT during one of the Her X-1 main-on states is 
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Fig. 3. Pulse period P of Her X-1 measured 
with Swift/BAT as a function of time. The 
cones around each point indicate the allowed 
range of the slope corresponding to the mea- 
sured P and its uncertainties. Measurement er- 
rors of the period itself are smaller than the 
symbol sizes. 



shown in Fig.|2l A straight Hne in the graph would correspond 
to a constant period defined by the slope of the line. The solid 
curve indicates the parabolic fit to the data corresponding to a 
constant positive P according to Eq. ^ (in the shown case the 
best-fit 1.4 X IQ-^^s/s). 



4. Pulse period variability 

Using the method described in the previous section we deter- 
mined the pulse period P and its time derivative P for most of the 
Her X-1 main-on cycles observed with BAT (for other parts of 
the 35 d cycle the flux was too low for such determinations). At 
the time of writing this paper the data are available for the time 
period from March 2005 to May 2009 that covers 45 35d cycles 
(main-on states) of the pulsar. For several cycles the BAT ob- 
servations have relatively poor statistics due to gaps in the data. 
For such cycles only P, but no P could be determined. For the 
cycle with the turn-on at MJD ~54339.68 even the pulse period 
could not be found due to a gap in the data. However, during 
this cycle Her X-1 wa s observed with the INTEGRAL satellite 
jStaubert et al.ll2009bl) . So, we included the value of P measured 
with INTEGRAL in our data set (since the INTEGRAL obser- 
vations caught only the end of the main-on, P could not be mea- 
sured). 

For the main-on starting at MJD ~53577.14 only the pulse 
period P could be measured with BAT (no P). But the main-on 
was simultane ously observed with the INTEGRAL and RXTE 
observatories (iKlochkov et a l. 2008; Stau bert et"aLll2009b ). The 
value of P measured during these pointings is consistent (within 
Icr uncertainty) with the value found from the BAT data. We 
consider this agreement as a successful cross-check between the 
Swift, RXTE, and INTEGRAL observations. 

The time evolution of the measured pulse period for the 
main-ons of Her X-1 is shown in Fig. [3] Where a corresponding 
value of P was measured, the 1 cr-uncertainty range is indicated 
by the cones, the orientation of which reproduce the measured P 
value. 

A close quantitative inspection of all measured values of P 
and P leads to the following list of statements about their evolu- 
tion with time. 

- The pulse period evolution shown in Fig.[3]resembles a saw- 
tooth where one can distinguish five spin-down and five spin- 
up episodes (from comparing adjacent measurements). The 



overall mean spin-up rate is < P >= -1.8203(3) x lO^'-'s s"' 
or -15.728(2)nsd"'. We note that this is slightly steeper 
than the general mean spin-up of -9nsd"' observed from 
the discovery of the source to the dramatic spin-down event 
during the Anomalous Low of 1999/2000 ( Klochkov 2007i) . 

- The spin-down episodes (with mean spin-down rates be- 
tween 0.95 X 10"'^ and 4.5 x 10"'^ s s"'), lasting from 1 to 
3 35 d cycles, are generally shorter than the spin-up episodes 
(with mean spin-up rates from 2.9 x lO"'-' to 4.4 x lO"'-' s s"') 
which last from 3 to 13 35 d cycles. 

- Locally measured P values, ranging from -1 1 x lO"'"' s s"' 
to +48 X lO^'^* s s"', show stronger modulation than P values 
found from comparing the pulse period of adjacent 35 d cy- 
cles which range from -6.0 x 10"'^ s s"^ to -h14x 10"^^ s s"^ 

- Generally, the locally measured P values and the "cycle-to- 
cycle" P values show a similar behavior (with the former 
having the larger amplitude). For many points, however, the 
locally measured P is substantially different from the one de- 
rived from P values of adjacent measurements (that is from 
the slope of the pulse period development). 

From the last two statements one can conclude that strong pulse 
period variations in Her X-1 occur on shorter time scales than 
the 35 d super-orbital period of the system. 

We will discuss the pulse period variations again in 
Section|6l 

5. Correlation between spin-up rate and X-ray flux 

A positive correlation between the spin-up rate of a neutron star 
and its X-ray l uminosity is generally expected from accretion 
theory (see e.g. lPringle & Reeslll972l) . where an increase of the 
mass accretion rate M leads to an increase in the rate of change 
of the neutron star's angular momentum. Such a correlation has 
been observed in many accreting pulsars (see e.g. lBildsten et al.l 
1 19971 for a review). In Her X-1, however, the correlation was 
so far questionable (see however IWilson et alj|1994l : iKlochkovl 
2007). One of the difficulties is the relatively low amplitude of 
the M variation in the system (reflected by the maximum main- 
on flux, see below), of a factor ~2. Another problem is the lack 
of local spin-up measurements (so far P in Her X-1 was locally 
measured only in a few cases). For the rest of the existing data 
P was estimated by taking differences of the pulse periods be- 
tween adjacent main-on states which gives an averaged value 
over two or more 35 d cycles (see e.g. IKlochkovl l2007l where 
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data from CGRO/BATSE were used). The X-ray luminosity Lx 
which we assume to be proportional to the mass accretion rate 
M is usually estimated from the maximum X-ray flux dur ing the 
main-on coincident with the observation (iStaubert et alJl2007h . 
Such an estimate is supposed to be very close to the local value 
of Lx- 

In this work we have used the values of P locally measured 
with BAT to explore the correlation be tween the spin-u p rate and 
the X-ray luminosity. Following Sta ubert et al.l (l2007h we have 
used the maximum main-on flux detected with BATQ as a mea- 
sure of Lx- In Fig.|4]we plot the measured pulse period derivative 
versus the maximum main-on flux. The dashed line indicates a 
linear fit to the data that takes uncertainties of both v ariables into 
account (using the orthogonal regression method, iBoggs et"aLl 
1198% . The data indicate an anticorrelation as predicted by the 
basic accretion theory. Inspection of the linear Pearson's corre- 
lation coefficient gives the probability of ~4xl0""^ to find the 
measured correlation by chance. We note, however, that such a 
high significance appears mainly due to the group of four points 
with high spin-down rate and low flux (in the upper left part of 
the graph in Fig.|4|i. The rest of the points form an uncorrected 
"cloud" around P - Q. On the other hand, the data contain no 
points with large spin-down [P > IQ^^s/s) and large luminos- 
ity (in the upper right part of the graph) that could destroy the 
correlation. 

The best-fit slope of the P(Lx) dependence assuming a lin- 
ear relation if found to be AP_i3/AFbat - -(1-2 ± 0.2) x 10^ 
where P^n = P/{10^^^ s/s) and Fbat is the count rate mea- 
sured with BAT. Assuming a distance of 7 kpc (see Introduction) 
one can find that 0.05 BAT cts s ' approximately corresponds to 
2x10-'^ erg s"' . Thus, we can rewrite the found value of the slope 
as AP^is/ALb? = -30 + 5 where L37 = Lx/(10")ergs-^ At a 
first approximation the value of the slope can be calculated using 
the equation: 



did) ■ I— u.^ 

— ^M,,^JGMRd-Ko^ + 

ut 
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1 



where the first term stands for the spin-up from the inner disk 
radius Rd, the second term describes spin-down due to magnetic 
coupling beyond the corotation radius Rc, and the last term takes 
into account angular momentum exchange due to magnetic cou- 
pling in the region between Rc and Rj, K() and k\ are numerical 
coefficients. The inner disk radius is about the Alfven radius of 
the mag netosphere ^a, wh ich is determined in the standard way 
(e.g. iPringle & Rees|[T972h : 



Ra = 



2M y/2GM 



2/7 



(4) 



At the corotational regime when Rj ~ R^ the last term in Eq. ^ 
can be rewritten in the form Kifi^3{Rc - Rd)/R'^ and is much 
smaller than the first two terms separately, but can be compara- 
ble to their difference at the point of equilibrium, where Rd - Rc- 
Indeed, from Eq. Q one can find that the spin-up and spin-down 
torques separately would give P — t- or - 10"" s s"' , whereas on 
average we observe < P >~ 10 s s The Equation. ^ can 
be rewritten as 

P^-^J^M^/GMR,-Ko^y (5) 



^ We used the Swift/BAT transient monitor results provided by the 

Swift/BAT teani 

( ,http://heas arc.gsfc.nasa.gov/docs/swift/results/transients/ 1 ). 
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Fig. 4. The locally measured time derivative of the pulse period, 
P, in Her X- 1 versus the maximum X-ray flux during the main- 
on state as determined from the Swift/BAT data. A linear fit (tak- 
ing uncertainties of both variables into account) is shown with 
the dashed line. 

and used to calculate the theoretical slope of the P{Lx) depen- 
dence. Here we did not include the third term in Eq. (O, which 
is by {Rc - Rd)/Rd times smaller than the spin-up term, and 
can be ignored when calculating the expected slope of P - Lx 
relation. Assuming a radiative efficiency of accretion of ~IQ% 
(Lx = O.lMc^) we can obtain dP-n/dL^i — -50 which is close 
to the value of AP^ij,/ AL^j found from the fit of the BAT data 
(-30 ± 5). A detailed discussion of the observed correlation will 
be given in Sects. l6^ and l6.3l 

6. Discussion 



(3) 6.1. Pulse period variations. Comparison witii BATSE data 



The general behaviour of the pulse period in Her X-1 measured 
with Swift/BAT (Fig.[3]l is similar to that observed previously 
in this source. The long-term spin-up trend is occasionally in- 
terrupted by short (a few 35 d cycles) spin-down episodes. In 
Fig.|5]we plot for comparison the historical development of the 
Her X-1 pulse period measured with CGRO/BATSE. We have 
made use of the publicly distributed p ulsar data as well as lists 
kindly provided by R. Wilson (see also lStaubert et al.ll2006[) . 

In both, the BATSE and BAT data, the spin-down episodes 
are generally shorter than the spin-up episodes. One also sees 
that in most cases the averaged absolute value of P is larger dur- 
ing the spin-down intervals than during spin-ups. Such an asym- 
metry is difficult to explain within the basic accretion theory 
(i.e. using Eq.|3]l assuming that the spin period reflects stochastic 
variations of M. In some cases spin-down episodes with subse- 
quent spin-ups are reminiscent of glitch-like behaviour observed 
in some radio-pulsars where a rapid change of the pulse pe- 
riod is followed by a slower "relaxation" to the long-term trend. 
However, compared to radio pulsars, the "glitches" in Her X-1 
appear with the opposite sign: an initial rapid increase of the pe- 
riod is followed by a slower decrease (as it should be if in both 
radio pulsars and Her X- 1 the glitches are caused by a sudden de- 
crease/increase of the moment of inertia of the star's crust). Two 
such glitch-like episodes are indicated in Fig.|5]by solid lines 
which show the fit of the "relaxation" intervals by an exponen- 
tial function on top of the linear decrease of the period (see e.g. 
IShemar & LvnelI996l) . The characteristic relaxation times, r, for 
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Fig. 5. Pulse period P of Her X-1 measured 
with CGRO/BATSE as a function of time. The 
values are derived from the publicly available 
pulse search data (R. Wilson, priv. comm.). The 
solid lines indicate the exponential fit to the 
pulse period during the "relaxation" time after 
two "glitch"-like episodes. 



the two cases are around 1 00 days which is simila r to the values 
observed in radio pulsars dShemar & Lvnel[l996l) . Some of the 
other spin-up/spin-down transitions also resemble glitches, even 
though their decay is less clearly exponential. One should note, 
however, than in radio pulsars the glitch-like variations of the 
spin period are predominantly observed in cases of young neu- 
tron stars, with characteristic ages (r = P/2P) less than 1 Myr. 
The neutron star in Her X- 1 is believed to be much older. 



6.2. Correlation P - Lx 

As we have shown, the BAT data indicate a correlation between 
the spin-up rate (-P) and the X-ray flux Lx with the slope of 
the coiTelation consistent with the prediction of the basic ac- 
cretion theory. However, even for the relatively narrow range of 
Lx observed in Her X-1 the correlation appears to be surpris- 
ingly loose. There are several physical and observational factors 
that might introduce dispersion in the observed P - Lx relation. 
The observational factors include the systematics that affects our 
measurements of the spin-up rate and the X-ray flux. As it was 
already mentioned in Sect.|5] the pattern-matching technique re- 
quires a stable shape of the pulse profiles. Even though we tried 
to restrict our analysis with the intervals where the pulse profile 
does not change significantly, we cannot exclude that the vari- 
ability of the profile shape contributes to the measured value of 
P. We believe, however, that this systematic effect, if present, 
does not exceed the statistical eiTors. On the other hand, the X- 
ray luminosity Lx might not be unambiguously related to the 
maximum main-on flux if the latter is subject to variable absorp- 
tion by a hot corona even in the middle of a main-on where the 
angle between the plane of the accretion disk and observer's line 
of sight reaches a maximum. 

A physical reason for the scattering in the correlation is a 
possible bimodal behaviour of the pulsar where it switches be- 
tween the spin-up and spin-down regimes at the same level of 
X-ray luminosity. Such pos sibility appears e . g. in the accretion 
torque theory presented by iLovelace et al.l (1 19951) where spin- 
up/spin-down transitions occur if the corotation radius is close 
to the Alfven radius (or, more precisely, to the turnover radius 
rto introduced by Lovelace et al. 1995), i. e. if the pulsar is close 
to the so-called corotational regime when the spin-up and spin- 
down torques are nearly balanced. 



The "uncorrected" cloud of data points around P = in 
Fig. |4] might indicate the area of the spin-up/spin-down transi- 
tions (bimodal behaviour) where different values of P are ob- 
served at the same level of Lx- 

Another factor that might increase the dispersion of data 
points in the correlation relates to the idea of a freely precessing 
neutron star in Her X-1 (Triimper et al. 1971). H aving a number 
of difficulties fr om the theoretical point of view (IShahamlll977l 
ISedrakian et al.|[T999L) . this hypothesis, however, can explain the 
syste matic variation of X-r ay pulse profiles exhibited by the sys- 
tem (IStaubert et al.ll2069 al and references therein). Free preces- 
sion results in systematic variations of the observed spin pe- 
riod of flie neutron star (IShakural ll988t iPostnov et all 119911; 
iBisnovatvi-Kogan & Kahabkal 19931) . With the geometrical pa- 
rameters of free precession that can be assumed for Her X-1 
on the basis of modelling of its pulse profiles (K. Postnov, priv. 
comm., paper is in preparation) the amplitude of P variation due 
to free precession might be as high as 10"'^ s/s which is compa- 
rable to the variations that we observe with BAT (Fig.|4|i. Since 
the precessional phase might be different in different observa- 
tions, the described effect might introduce additional scattering 
in the P-Lx correlation if the free precession indeed takes place 
in Her X-1. 



6.3. Evidence for ttie coronal mass ejection 

As mentioned in Sect.|5] in Her X-1 spin-up and spin-down 
torques are very well balanced, so that the averaged P is rela- 
tively close to zero. However, from the points in the upper left 
part of the graph in Fig.|4]one can conclude that there occur ex- 
tremely large spin-down torques (at small fluxes), which are up 
to 5 times as strong as spin-up ones. So, the key feature to ex- 
plain is why do we observe such strong spin-down episodes? 

One possibility is to assume that the accretion disk carries 
some magnetic field which can interact via reconnection with the 
neutron star's magnetosphere beyond the corotation radius Rc. 
This might imply that beyond Rc the field lines can sometimes 
inflate to become open (see also discussion in ILovelace et al.l 
1995). During such episodes, a substantial fraction of matter in 
the inner part of the accretion disk can escape the system in the 
form of a coronal wind ejection along the open field lines. 
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Such an ejection of matter should be reflected in a secu- 
lar change of the system's orbital period which is ind eed ob - 
served in Her X-1 (Deeteretal. 1991; Staubert et al. 2 009bh . 
To assess the importance of coronal mass ejections for the or- 
bital period evolution, we invoke general considerations of the 
non-conservative treat ment of binary orbital parameters (see e.g. 
iGrishchuk et al.ll200lh . We shall assume a circular binary orbit. 
Let Mo be the mass of the optical star and - that of the neu- 
tron star, q = M^/Mg is the binary mass ratio. The total angular 
momentum of the binary is mostly stored in the orbital motion 
of the binary components: J = /orb = M^-Mg/iMf + M„)(x)orbfl^, 
where Worb = ^n/Porb is the orbital angular frequency, Poit is the 
orbital period, and a - + is the binary's semimajor axis. 
Let us de fine the non-conserva tiveness parameter in the standard 
way (e.g. lRitter&Kolblll992l) : 



77= r-^ < 1, Mo<0. 



(6) 



Next we assume that the ejected mass carries away the spe- 
cific angular momentum of the neutron star jx - Woibfl^ - 
(jj oihO^ (Mo I {Mx + M„yf-. (This assumption is justified by the 
small inner radius of the disk <sc a.v). Then using the Kepler's 
3-rd law and the total angular momentum balance, J - jAMx + 
Mo), we can express the fractional change of the orbital period 
through MxIMx, q and 77: 



I f orb 
3Porb 



1-^ 



^ l\ g/3 + l 
T]] q+\ 



(7) 



For ?7 = 1 the familiar expression for the conservative mass ex- 
change is obtained. It is convenient to normalize the mass accre- 
tion rate onto the neutron star M, to the value that can be derived 
from the observed fractional change of the orbital period in the 
conservative case: 



1 



M, 



3 ^"01* - 1 ' 



(8) 



s. In Her X-1 q ^ 0.63 and one 



8 X 10' ^ gs ' for th e measured /"orb 



so that m s M.,/(M,)cons- 

finds (M,)cons ., 

-4.85 X 10 " ss ' (IStaubertet al.ll2009bb . Then we can elim- 
inate PoTb/Porb from the left hand side of Eq. (jT) to obtain the 
equation for 77 at a given m: 



q^ + \q-\ 



(9) 



From here we see that m < 1 leads to 77 < 1, i. e. if one wants 
to decrease the mass accretion rate onto the neutron star to get 
smaller X-ray luminosity (as is the case of Her X-1, where the 
mean observed ~ 2 X 10^^ ergs"' is 4 times smaller than 
the one following from t he conservative mass exchange analy- 
sis, IStaubert et al.ll2009 b). a non-conservative mass exchange is 
required. Specifically, if we want to bring in accordance the ob- 
served Poib and Lx in Her X-1, we would need m - 1/4 and 
(from Eq.|9]l 77 ~ 0.1, a fairly high non-conservative mass ex- 
change. 

Here we should note that the binary mass ratio q in Her X-1 
may be uncertain. Indeed, recent analysis of non-LTE effects in 
the formation of the Hy absorption line allow s two solutions: 
q = 0.45 and q = 0.72 (Abubekerov et al.ll2008l) . Eq. © implies 
that the critical (M^)cons decreases for smaller q: for q - 0.45 
we obtain (MiOcons - 5.4 x 10'^ gs"'. For this mass ratio the 
appropriate value of the dimensionless parameter m is 1/2, for 



which from Eq. (|9j we find 77 ^ 0.38. This value implies that 
on average about half of the matter transferred through the disk 
should escape from the system to provide the observed decrease 
rate of the system's orbital period. Note, however, that such a 
small mass ratio suggests an unusually small mass of the neutron 
star about 0.85 Mq. 

It is very likely that an accretion disk corona is present in 
Her X- 1 (see recent analysis of Chandra X-ray observations 
iJi et al.| [2009). This suggests that there might be a permanent 
coronal accretion disk wind which carries out some angular mo- 
mentum from the system. Episodic ejection of matter in Her X-1 
has sporadically app eared in the literature to e xplain some prop- 
erties of the svstem. lCrosa & Bovntonl(ll980l) found that the av- 
eraged mass transfer rate to the outer rim of the accretion disk is 
somewhat larger than that required to maintain the obs erved X- 
ray luminosity. In the model of Scha ndl & Meverld 19941) the disk 
wind results from irradiatio n of the dis k by the central source. 
IVrtilek et al] (12001) and Boroson et al.1 dMol,) found signatures 
of outflowing gas in the UV spectrum of the system. In the frame 
of our model, the mass ejection from the system through the 
open magnetic field lines occurs most efficiently during strong 
spin-down episodes which are associated with small X-ray lu- 
minosity. Indeed, as we see in Fig|4l the observed X-ray flux 
is decreased by a factor of two during strong spin-down. From 
Eq. (|9]l it is easy to find that at a given q a fractional decrease 
in m leads to comparable fractional decrease in 77, i.e. accretion 
indeed becomes more non-conservative during strong mass ejec- 
tion episodes. During such episodes, the neutron star spin-down 
power luio is spent to expel accreting matter from the inner disk 
radius R^ ~ Rc'. 

GM 



Icocj ~ M, 



Rc 



(10) 



This equation is satisfied for the observed parameters of 
Her X-1: ejected mass rate during strong spin-downs Mej ~ 
0.5Mi. ^ gs-\P^ lO-'^ss"', and/?,. - 1.3 x 10** cm. 

7. Summary 

We have used the publicly available Swift/BAT data to study 
the long-term behaviour of the pulsar's spin period in Her X-1. 
The measured pulse period variations were compared with those 
observed previously with CGRO/BATSE. 

For the first time the pulse period derivative was measured 
for a long series of observed main-on states of the source. This 
allowed us for the first time to test the correlation between the 
X-ray luminosity and the locally measured spin-up rate of the 
neutron star in Her X- 1 . We argue that the data indeed show such 
a correlation with the slope consistent with the prediction of the 
basic accretion theory for the parameters of Her X-1. The rela- 
tively large scattering of the data point in the vicinity of P ~ 
can be caused by the bimodal behaviour of the accretion flow 
configuration at the magnetospheric boundary which results in 
switching of the pulsar between spin-up and spin-down branches 
when it stays close to the corotational regime. In addition, free 
precession of the neutron star, if it takes place in the system, 
can affect significantly the measured values of P and, therefore, 
contribute to the scattering. 

We argue that together with the long-term decrease of the 
orbital period in Her X-1 the measured pulse period behaviour 
requires the presence of mass ejection from the inner parts of 
the accretion disk along the open magnetic field lines. The mass 
ejection episodes probably take place during strong spin-down 
episodes which are associated with the small X-ray luminosity. 
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We would like to note that the described technique al- 
lows one to use the BAT instrument as a long-term moni- 
tor of spin periods in othe r bright accreting pulsars (see also 
ICamero-Arranz et al.ll20()9h . 
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